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ABSTRACT 

We  report  on  the  physical  properties  of  solar  sequential  chromospheric  brightenings  (SCBs)  observed  in  conjunction 
with  moderate- sized  chromospheric  flares  with  associated  Coronal  mass  ejections.  To  characterize  these  ephemeral 
events,  we  developed  automated  procedures  to  identify  and  track  subsections  (kernels)  of  solar  flares  and  associated 
SCBs  using  high-resolution  Hck;  images.  Following  the  algorithmic  identification  and  a  statistical  analysis,  we 
compare  and  find  the  following:  SCBs  are  distinctly  different  from  flare  kernels  in  their  temporal  characteristics 
of  intensity,  Doppler  structure,  duration,  and  location  properties.  We  demonstrate  that  flare  ribbons  are  themselves 
made  up  of  subsections  exhibiting  differing  characteristics.  Flare  kernels  are  measured  to  have  a  mean  propagation 
speed  of  0.2  km  s“^  and  a  maximum  speed  of  2.3  km  s“^  over  a  mean  distance  of  5  x  10^  km.  Within  the  studied 
population  of  SCBs,  different  classes  of  characteristics  are  observed  with  coincident  negative,  positive,  or  both 
negative  and  positive  Doppler  shifts  of  a  few  km  s“^  The  appearance  of  SCBs  precedes  peak  flare  intensity  by 
~12  minutes  and  decay  ~1  hr  later.  They  are  also  found  to  propagate  laterally  away  from  flare  center  in  clusters 
at  45  km  s“^  or  117  km  s“^  Given  SCBs’  distinctive  nature  compared  to  flares,  we  suggest  a  different  physical 
mechanism  relating  to  their  origin  than  the  associated  flare.  We  present  a  heuristic  model  of  the  origin  of  SCBs. 

Key  words:  Sun:  chromosphere  -  Sun:  flares 


1.  INTRODUCTION 

The  solar  chromosphere  exhibits  three  different  classes  of 
small-scale  intensity  brightenings:  flare,  plage,  and  compact 
brightenings.  Although  each  is  characterized  by  an  enhanced 
temporal  Ha  brightness  relative  to  a  background  quiet  Sun,  they 
each  have  distinct  physical  processes  governing  their  spatial  and 
temporal  evolution.  Typically,  brightenings  have  been  identified 
and  characterized  manually  from  a  single  data  source  (e.g.,  Kurt 
et  al.  2000;  Ruzdjak  et  al.  1989;  Veronig  et  al.  2002).  However,  in 
order  to  form  a  better  understanding  of  the  underlying  dynamics, 
data  from  multiple  sources  must  be  utilized  and  numerous 
similar  features  must  be  statistically  analyzed. 

This  work  focuses  on  flare  brightenings  and  associated  com¬ 
pact  brightenings  called  sequential  chromospheric  brightenings 
(SCBs).  SCBs  were  first  observed  in  2005  and  appear  as  a  series 
of  spatially  separated  points  that  brighten  in  sequence  (Bala- 
subramaniam  et  al.  2005).  SCBs  are  observed  as  multiple  trains 
of  brightenings  in  association  with  a  large-scale  eruption  in  the 
chromosphere  or  corona  and  are  interpreted  as  progressive  prop¬ 
agating  disturbances.  The  loci  of  brightenings  emerge  predom¬ 
inantly  along  the  axis  of  the  flare  ribbons.  SCBs  are  correlated 
with  the  dynamics  which  cause  solar  flares,  coronal  restructur¬ 
ing  of  magnetic  fields,  halo  Coronal  mass  ejections  (CMEs), 
EIT  waves,  and  chromospheric  sympathetic  flaring  (Balasub- 
ramaniam  et  al.  2005).  Pevtsov  et  al.  (2007)  demonstrate  that 
SCBs  have  properties  consistent  with  aspects  of  chromospheric 
evaporation. 

This  article  presents  a  new  description  of  the  dynamical  prop¬ 
erties  of  SCBs  resulting  from  applying  a  new  automated  method 
(Kirk  et  al.  2011)  of  identifying  and  tracking  SCBs  and  associ¬ 
ated  flare  ribbons.  This  tracking  technique  differs  from  previous 
flare  tracking  algorithms  in  that  it  identifies  and  tracks  spatial 
and  temporal  subsections  of  the  flare  and  all  related  brightenings 
from  pre-flare  stage,  through  the  impulsive  brightening  stage. 


and  into  their  decay.  Such  an  automated  measurement  allows  for 
tracking  dynamical  changes  in  intensity,  position,  and  derived 
Doppler  velocities  of  each  individual  subsection.  The  tracking 
algorithm  is  also  adapted  to  follow  the  temporal  evolution  of 
ephemeral  SCBs  that  appear  with  the  flare.  In  Section  2,  we  de¬ 
scribe  the  data  used  to  train  the  algorithm  and  the  image  process¬ 
ing  involved  in  the  detection  routine.  In  Section  3,  we  present  the 
results  of  tracking  the  evolution  of  flare  kernels  through  an  erupt¬ 
ing  flare.  In  Section  4,  we  present  the  application  of  the  tracking 
algorithm  to  ephemeral  SCBs.  We  present  a  physical  interpreta¬ 
tion  of  the  SCBs  and  provide  a  heuristic  model  of  the  origin  of 
SCBs  in  Section  5.  Finally,  in  Section  6  we  discuss  implications 
of  these  results  and  provide  future  direction  for  this  work. 

2.  DATA  AND  DATA  PROCESSING 

In  this  study,  we  use  chromospheric  Hck;  (6562.8  A)  images 
from  the  USAF’s  Improved  Solar  Observing  Optical  Network 
(ISOON;  Neidig  et  al.  1998)  prototype  telescope  to  study  flare 
ribbons  and  SCBs.  ISOON  is  an  automated  telescope  producing 
2048x2048  pixel  full-disk  images  at  a  one-minute  cadence. 
Each  image  has  a  V'.  1  spatial  sampling,  is  normalized  to  the  quiet 
Sun,  and  corrected  for  atmospheric  refraction  (Figures  1  (a),  2(a), 
and  3(a)).  Nearly  coincident  to  the  spectral  line  center  images, 
ISOON  also  records  Ha  zb  0.4  A  off-band  Doppler  images,  from 
which  Doppler  signals  are  derived. 

For  this  study,  we  chose  to  apply  the  brightening  detection 
algorithms  to  three  flares  where  Balasubramaniam  et  al.  (2006) 
had  previously  identified  SCB  events  (Table  1).  Each  of  the 
events  selected  has  a  two-ribbon  configuration  and  an  associated 
halo  CME.  Of  these,  both  the  2005  May  6  and  2005  May 
13  events  were  both  located  near  disk  center,  while  the  2004 
November  9  event  was  near  the  western  limb.  Images  were 
extracted  from  the  archive  from  ±3.5  hr  of  the  eruption  start 
time,  yielding  a  data  cube  with  ~400  images  for  each  event. 
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Figure  1.  Two-ribbon  flaring  event  from  2005  May  13.  Panel  (a)  shows  an  example  of  a  calibrated  Ho;  ISOON  image  with  the  ROI  highlighted.  The  ISOON  images 
are  2048  x  2048  pixels  with  each  pixel  having  a  T.'l  resolution.  Panel  (b)  is  the  ROI  after  preprocessing.  This  ROI  covers  576  x  576  pixels  corresponding  to 
2.12  X  10^^  km^  on  the  solar  surface.  Panel  (c)  is  a  Doppler  measure  of  the  same  ROI  as  part  (b).  The  Doppler  velocity  image  is  created  out  of  ISOON  Ho;  off-band 
images  described  in  Section  2  where  the  velocity  ranges  from  —26.6  to  21.5  km  s“^  from  black  to  white.  Panel  (d)  shows  intensity  curves  over  the  time  period  of 
interest:  13:00-20:02  UT  as  described  in  detail  in  the  text. 


Table  1 

Events  Selected  for  Study 


Date 

Start  Time 
(UT) 

Duration 

(hr) 

Flare  Class 

CME 

2004  Nov  9 

16:59 

0.5 

M8.9 

Halo 

2005  May  6 

16:03 

2.1 

C8.5 

Halo 

2005  May  13 

16:13 

1.3 

M8.0 

Halo 

Note.  Each  event  was  visually  identified  as  having  a  two-ribbon  configuration 
and  SCBs  associated  with  the  flare. 


An  Ha  relative  intensity  curve,  an  Ha  maximum  intensity 
curve,  and  Geostationary  Operational  Environmental  Satellites 
{GOES)  hard  and  soft  X-ray  fluxes  are  plotted  to  characterize 
the  flare  (Figures  1(d),  2(d),  and  3(d))  as  described  in  Kirk  et  al. 
(2011). 

Each  ISOON  image  is  reconditioned  to  remove  solar  limb 
darkening.  The  images  are  then  de-projected  into  conformal 
coordinates  using  a  Guyou  projection  (an  oblique  aspect  of  the 
Peirce  projection;  Peirce  1879),  which  removes  the  projection 
effects  of  imaging  the  solar  sphere.  Each  image  is  then  cropped 
to  the  region  of  interest  (ROI;  Figures  1(b),  2(b),  and  3(b)).  The 
set  of  images  are  aligned  using  a  cross-correlation  algorithm 
eliminating  the  rotation  effects  of  the  Sun.  Frames  that  contain 
bad  pixels  or  excessive  cloud  interference  are  removed. 

The  red  and  blue  wings  of  ISOON  Doppler  images  are  each 
preprocessed  using  the  same  technique  as  applied  to  the  line 


center  images  described  above  (Figures  1(c),  2(c),  and  3(c)).  In 
order  to  produce  velocity  measurements,  a  Doppler  cancelation 
technique  is  employed  in  which  red  images  are  subtracted 
from  blue  images  (for  a  modern  example,  see  Connes  1985). 
Typically,  red  and  blue  images  are  separated  by  about  4  s  and  no 
more  than  5  s.  A  mean  zero  redshift  in  the  subtracted  Doppler 
image  in  the  quiet  Sun  serves  as  a  reference.  In  dynamical 
situations  such  as  solar  flares,  Ha  profiles  are  often  asymmetric 
especially  in  emission  making  this  Doppler  technique  invalid. 
To  avoid  asymmetric  profiles,  only  values  derived  outside 
of  the  flaring  region  are  considered  such  that  the  areas  of  interest 
at  ±0.4  A  are  still  in  absorption  even  if  raised  in  intensity 
(Balasubramaniam  et  al.  2004).  The  values  from  the  subtracted 
images  are  translated  into  units  of  km  s“^  by  comparison  to  a 
measured  response  in  spectral  line  shifts,  which  are  calibrated 
against  an  intensity  difference  for  the  ISOON  telescope.  In  this 
context,  we  interpret  the  Doppler  shift  as  a  line-of- sight  velocity. 
Thus,  Doppler  velocity  [f/)]  is  defined  as 

—  k{^hhjiQ  ^^red)?  (1) 

where  I  is  the  measured  intensity  and  ^  is  a  linear  fitted 
factor  that  assumes  the  intensity  changes  due  to  shifts  in  the 
symmetric  spectral  line,  which  can  be  attributed  to  a  Doppler 
measure  as  a  first-order  approximation.  The  linear  factor  [k] 
has  been  independently  determined  for  the  ISOON  telescope 
by  measuring  the  full  Ha  spectral  line  Doppler  shift  across 
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Figure  2.  Two-ribbon  flaring  event  from  2005  May  6.  Panel  (a)  shows  an  example  of  a  calibrated  Ho:  ISOON  image  with  the  region  of  interest  (ROI)  highlighted. 
The  ISOON  images  are  2048  x  2048  pixels  with  each  pixel  having  a  17 1  resolution.  Panel  (b)  is  the  ROI  after  preprocessing.  This  ROI  coves  768  x  768  pixels 
corresponding  to  3.77  x  10^^  km^  on  the  solar  surface.  Panel  (c)  is  a  Doppler  measure  of  the  same  ROI  as  part  (b).  The  Doppler  velocity  image  is  created  out  of 
ISOON  Ho:  off-band  images  described  in  Section  2  where  the  velocity  ranges  from  —12.5  to  6.9  km  s“^  from  black  to  white.  Panel  (d)  shows  intensity  curves  over 
the  time  period  of  interest:  12:28-22:47  UT  as  described  in  detail  in  the  text. 


the  entire  solar  disk  (K.  S.  Balasubramaniam  2011,  private 
communication) . 

2.1.  Detection  and  Tracking 

An  animated  time  series  of  sequential  images  covering 
an  erupting  flare  reveals  several  physical  characteristics  of 
evolving  ribbons:  the  ribbons  separate,  brighten,  and  change 
their  morphology.  Adjacent  to  the  eruption,  SCBs  can  be 
observed  brightening  and  dimming  in  the  vicinity  of  the  ribbons. 
Kirk  et  al.  (201 1)  describe  in  detail  techniques  and  methods  used 
to  extract  quantities  of  interest  such  as  location,  velocity,  and 
intensity  of  flare  ribbons  and  SCBs.  The  thresholding,  image 
enhancement,  and  feature  identification  are  tuned  to  the  ISOON 
data.  The  detection  and  tracking  algorithms  are  specialized  for 
each  feature  of  interest  and  require  physical  knowledge  (e.g., 
size,  peak  intensity,  and  longevity)  of  that  feature  being  detected 
to  isolate  the  substructure. 

Briefly,  the  detection  algorithm  first  identifies  candidate 
bright  kernels  in  a  set  of  images.  In  this  context,  we  define  a 
kernel  to  be  a  locus  of  pixels  that  are  associated  with  each  other 
through  increased  intensity  compared  with  the  immediately 
surrounding  pixels.  Each  kernel  has  a  local  maximum,  must 
be  separated  from  another  kernel  by  at  least  one  pixel,  and  does 
not  have  any  predetermined  size  or  shape.  Next,  the  algorithm 
links  time-resolved  kernels  between  frames  into  trajectories.  A 
filter  is  applied  to  eliminate  inconsistent  or  otherwise  peculiar 


detections.  Finally,  characteristics  of  bright  kernels  are  extracted 
by  overlaying  the  trajectories  over  complementary  data  sets.  To 
aid  in  this  detection,  tracking  software  developed  by  Crocker, 
Grier,  and  Weeks  was  used  as  a  foundation  and  modified  to  fit 
the  needs  of  this  project^  (Crocker  &  Grier  1996). 

In  order  to  characterize  a  kernel,  we  calculate  its  integrated 
intensity,  radius  of  gyration,  and  eccentricity.  The  eccentricity 
of  the  kernel  (as  defined  by  its  semi-major  {a)  and  semi-minor 
(h)  axis)  is  calculated  using 


The  integrated  intensity  of  a  given  kernel  is  defined  by 

m  =  ^  A{x  +  i,y  +  j),  (3) 

where  A(x  i,  y  j)  is  the  intensity  of  the  pixel  located  at 
(x-\-i,y-\-  j),  the  kernel  has  a  brightness-weighted  centroid  with 
coordinates  (x,  y),  and  co  is  the  radius  of  the  mask  (Crocker  & 
Grier  1996).  The  radius  of  the  mask  is  chosen  to  be  8  pixels  in 
this  case.  A  radius  of  gyration,  Rg,  is  related  to  the  moment  of 


^  Crocker’s  software  is  available  online  at  www.physics.emory.edu/ 
~weeks/idl/. 
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Figure  3.  Two-ribbon  flaring  event  from  2004  November  9.  Panel  (a)  shows  an  example  of  a  calibrated  Ho;  ISOON  image  with  the  ROI  highlighted.  The  ISOON 
images  are  2048  x  2048 pixels  with  each  pixel  having  a  l'!l  resolution.  Panel  (b)  is  the  ROI  after  preprocessing.  This  ROI  covers  600  x  600 pixels  corresponding  to 
2.30  X  10^^  km^  on  the  solar  surface.  Panel  (c)  is  a  Doppler  measure  of  the  same  ROI  as  part  (b).  The  Doppler  velocity  image  is  created  out  of  ISOON  Ho;  off-band 
images  described  in  Section  2  where  the  velocity  ranges  from  —9.3  to  14.6  km  s“^  from  black  to  white.  Panel  (d)  shows  intensity  curves  over  the  time  period  of 
interest:  14:37-23:58  UT  as  described  in  detail  in  the  text. 


inertia  [/]  by  using 


I  =  '^mkrl  =  niRl,  (4) 

k 

r  =  {a  +  b)/2,  (5) 

where  is  the  mass  of  particle  k,  m  is  the  total  mass  of  the 
system,  and  r  is  the  distance  to  the  rotation  axis.  In  this  context, 
we  interpret  the  radius  of  gyration  to  be 


R 


2 

8 


(6) 


where  m  is  the  integrated  intensity  as  defined  in  Equation  (3) 
(Crocker  &  Grier  1996).  Figure  4  diagrams  a  hypothetical 
kernel.  The  results  of  this  algorithm  identify  and  characterize 
several  kernels  in  the  flaring  region  (Figure  5).  The  total  number 
of  flare  kernels  detected  is  typically  between  100  and  200,  while 
the  number  of  SCB  kernels  detected  is  typically  three  to  four 
times  that  number. 


3.  FFARE  RIBBON  PROPERTIES 

The  properties  of  kernels  identified  in  flare  ribbons  can  be 
examined  in  two  ways:  (1)  each  kernel  can  be  considered  as  an 
independent  aggregation  of  compact  brightenings  or  (2)  kernels 
can  be  considered  as  dependent  on  each  other  as  fragments 


of  a  dynamic  system.  Each  category  brings  about  contrasting 
properties.  If  one  considers  the  kernels  as  independent  elements, 
this  provides  a  way  to  examine  changes  to  subsections  of  the 
flare  and  is  discussed  in  Section  3.1.  Associating  kernels  with 
their  contextual  surroundings  allows  a  way  to  examine  the  total 
evolution  of  the  flare  without  concern  of  how  individual  kernels 
behave.  This  type  of  examination  is  addressed  in  Section  3.2. 

3.1.  Qualities  of  Individual  Kernels 

A  sample  of  six  kernels  from  the  2005  May  13  event,  letters 
U-X  in  Figure  6,  is  representative  of  the  majority  of  flare  kernels 
tracked.  The  kernels  were  selected  from  three  different  regions 
of  this  flare.  The  normalized  Ha  intensity  of  each  of  these 
kernels  is  shown  in  Figure  7.  In  contrast  with  the  integrated  flare 
light  curve,  an  individual  flare  kernel  often  has  a  shorter  lifetime. 
This  is  because  an  individual  kernel  may  first  appear  in  the 
impulsive  phase  of  the  flare  (as  demonstrated  in  kernels  W-Y), 
while  other  kernels  may  disappear  as  the  flare  begins  to  decay  (as 
in  kernels  U  and  Z).  Most  likely  these  kernels  merged  with  one 
another,  at  which  point  their  unique  identity  was  lost.  Kernels  W 
and  X  are  the  only  two  that  show  the  characteristic  exponential 
dimming  found  in  the  reference  curves  in  Figure  1(d).  All  of 
these  kernels  have  peak  intensities  within  a  few  minutes  of  the 
peak  of  the  total  flare  intensity  and  have  sustained  brightening 
above  background  levels  for  up  to  an  hour. 

The  integrated  speed  of  displaced  kernels  provides  context 
to  the  evolution  of  intensities.  We  define  the  integrated  speed 
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Figure  4.  Diagram  of  a  hypothetical  kernel  (either  flare  or  SCB)  used  for 
classiflcations  in  this  paper.  Each  pixel  is  labeled  with  a  representative  intensity 
[Ai].  The  semi-major  axis  [a]  and  semi-minor  axis  [b]  are  shown.  This  kernel 
has  the  following  properties:  total  intensity,  mo  =  43,  semi-major  axis,  a  =  4, 
semi-minor  axis,  b  =  2,  eccentricity,  e  =  0.87,  centroid,  (x,  y)  =  (0.26,  0.18), 
and  radius  of  gyration,  Rg  =  1.9.  The  maximum  intensity,  in  this  case  5,  is 
deflned  as  having  pixel  coordinates  (0,  0). 


of  displacement  to  be  the  sum  over  all  time  steps  of  the 
measured  velocity  of  an  individual  kernel  and  thus  small  velocity 
perturbations  are  minimized  by  the  sum.  The  peak  integrated 
speed  measured  for  each  kernel  peaks  at  ~2.3  km  s“^  and  has  a 
mean  of  ~0.2  km  s“^ .  There  is  significant  motion  along  the  flare 
ribbons  as  well  as  outflow  away  from  flare  center.  The  motions 
are  complex  but  generally  diverge.  The  total  spatial  (angular) 
displacement  of  each  flare  kernel  peaks  at  ~20.7  Mm  and  has 
a  mean  of  "^5.0  Mm.  Flare  kernels  that  exhibited  the  greatest 
speed  did  not  necessarily  have  the  greatest  displacements.  As 


Figure  6.  Image  from  the  2005  May  13  event.  Letters  A-F  point  to  the  locations 
several  SCB  kernels  that  are  shown  in  more  detail  in  Figure  7.  Letters  U-Z 
point  to  the  locations  of  a  sample  of  several  flare  kernels  that  are  also  shown 
in  more  detail  in  Figure  7.  Note  that  since  flare  kernels  move  with  the  evolving 
flare  ribbon,  the  locations  of  these  kernels  are  only  accurate  for  this  image.  In 
contrast,  SCB  kernels  are  stationary  and  remain  in  the  same  location  in  any 
image. 


the  flare  develops,  most  of  the  motion  in  flare  ribbons  indicated 
is  synchronous  with  the  separation  of  the  ribbons.  A  few  tracked 
kernels  indicate  motion  along  flare  ribbons.  These  flows  along 
flare  ribbons  are  consistently  observed  in  all  three  flares. 

3.2.  Derived  Flare  Quantities 

Examining  trajectories  of  flare  kernels  offers  insight  into 
the  motions  of  the  flare  ribbons  as  they  evolve  through  the 
eruption.  Figure  8  shows  a  time  series  of  images  centered 
about  the  peak  time  of  the  flare.  The  velocities  of  the  detected 
flare  kernels  are  superposed  as  vectors.  Note  the  initial  outflow 
of  the  two  ribbons  near  the  flare  peak  at  16:49  UT.  But,  as 
the  flare  continues  to  evolve,  there  is  more  motion  along  the 
flare  ribbons.  This  is  probably  a  result  of  the  overarching 
loops  readjusting  after  the  reconnection  event.  Beneath  each 
image  is  a  histogram  of  the  distribution  of  velocities  at  each 
time  step.  The  mean  apparent  lateral  velocity  of  each  kernel 
remains  below  0.5  km  s“^  throughout  the  flare  even  though 


Figure  5.  Example  ISOON  Ho:  images  from  the  2005  May  13  two-ribbon  flare.  The  detection  algorithm  is  applied  to  the  image  series  to  extract  either  ribbon 
brightening  or  compact  brightening.  The  left  image  shows  one  time  slice  of  individual  SCBs  identifled  external  to  the  flare  ribbon.  The  right  image  shows  the  same 
time  slice  with  brightenings  identifled  along  the  flare  ribbon.  Each  detection  is  made  up  of  a  small  matrix  of  pixels  associated  with  the  brightening  called  a  kernel. 
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Time  (UT) 


Figure  7.  Ho;  normalized  intensity  curves  for  each  of  the  kernels  indicated  in  Figure  6  by  the  associated  letters.  On  the  left,  the  tracked  flare  kernels  are  shown.  The 
normalized  intensity  curves  for  the  SCB  kernels,  shown  on  the  right,  are  stationary  and  therefore  tracked  throughout  the  flaring  window.  In  both  panels,  the  vertical 
dashed  line  indicates  the  time  at  which  the  flare  intensity  peaks  (16:49  UT)  derived  from  integrating  all  flare  kernel  intensities  at  each  time  step. 


16:36  16:46  16:56  17:07  17:17  17:27 
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Figure  8.  Time  series  of  images  with  the  apparent  lateral  velocity  of  each  tracked  flare  kernel  at  that  time  step  plotted  over  the  Ho;  flare  ribbons.  A  histogram  of  the 
number  of  flare  kernels  as  a  function  of  velocity  corresponding  to  each  time  step  is  plotted  underneath. 


the  total  velocity  changes  significantly.  Qualitatively,  the  bulk 
of  the  apparent  motion  appears  after  the  flare  peak  intensity 
at  16:49  UT. 

Grouping  the  kernels’  velocities  into  appropriate  bins  pro¬ 
vides  another  way  to  analyze  the  dynamics  of  the  flare  (Figure  9). 
The  kernels  with  apparent  lateral  speeds  above  0.4  km  s“^  are 
highly  coincident  with  the  peak  intensity  of  the  flare.  The  num¬ 
ber  of  these  fast-moving  kernels  peaks  within  a  couple  minutes 
of  flare  peak,  and  then  quickly  decays  back  to  quiet  levels.  The 
velocity  bins  between  0.1  and  0.4  km  s“^  show  a  different  sub¬ 
structure.  These  velocity  bins  peak  ~30  minutes  after  the  peak 
of  the  flare  intensity.  They  show  a  much  slower  decay  rate. 


staying  above  the  pre-flare  velocity  measured  throughout  the 
decay  phase  of  the  Hck:  flare  curve. 

Figure  10  (left  column)  shows  integrated  unsigned  kernel 
velocity  (speed)  for  each  image.  In  general,  the  evolution  of  the 
integrated  kernel  speed  has  a  shape  similar  to  the  intensity  curve. 
For  the  2005  May  6  event,  the  speed  curve  is  more  similar  to  the 
GOES  X-ray  intensity  curve  than  the  Hey  curve.  For  the  2004 
November  event,  the  limb  geometry  of  the  flare,  as  well  as  the 
noise  in  the  original  data  set,  results  in  a  noisy  curve.  Despite 
these  difficulties,  a  clear  increase  in  total  speed  is  apparent  near 
the  peak  of  the  flare.  A  continuum  integrated  speed  level  of 
1-2  km  s“^  is  consistent  in  each  event,  indicating  the  quality  of 
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Figure  9.  Number  of  elements  in  each  apparent  lateral  velocity  bin  plotted 
as  a  function  of  time.  The  vertical  dashed  line  marks  the  time  of  the  peak 
flare  intensity.  The  velocity  bins  are  0.1  km  s“^  in  width.  The  smallest  bin, 
0-0.1  km  s“^  is  excluded  because  of  a  base  level  of  noise  inherent  in  the 
detection  routine. 

data  and  of  the  tracking  software.  The  peak  integrated  speed  is 
about  15  km  s“^  for  all  three  flares. 

Examining  the  change  in  the  speed  of  each  kernel  be¬ 
tween  successive  images  yields  a  derived  acceleration  for  each 
tracked  kernel.  The  integrated  unsigned  acceleration  is  plotted  in 
Figure  10  (middle  column).  Both  the  2005  May  13  and  2005 
May  6  events  show  peak  acceleration  after  the  peak  of  the 
intensity  curve.  Hence,  the  majority  of  the  acceleration  in  the 
apparent  motion  of  the  flare  comes  in  the  decay  phase  of  Hck: 
emission,  in  concert  with  the  formation  of  post-peak  flare  loops, 
as  seen  in  movies  of  Transition  Region  and  Coronal  Explorer 
images  of  other  flares.  The  peak  acceleration  appears  uncorre¬ 
lated  with  the  strength  of  the  flare  since  the  peak  of  the  M  class 
May  13  event  is  just  over  300  km  s“^,  while  the  peak  of  the 
C-class  flare  on  May  6  is  nearly  400  km  s“^.  The  acceleration 
curve  for  the  November  9  flare  has  too  much  noise  associated 
with  it  to  decisively  determine  the  peak  value. 

The  right  column  of  Figure  10  shows  a  derived  the  kinetic 
energy  associated  with  the  measured  motion  of  the  flare  kernels. 
This  was  accomplished  by  assuming  a  chromospheric  density 
[a]  of  10“^  kg  m“^  and  a  depth  [h]  of  1000  km.  The  mass 
contained  in  each  volume  is  then  calculated  by 

m  =  iihoR^^\J  1  —  e^,  (7) 

where  is  defined  in  Equation  (6)  and  e  is  defined  in 
Equation  (2).  The  kinetic  energy  under  a  kernel  is  defined  in 
the  standard  way  to  be 

(8) 

where  v  (speed)  is  a  measured  quantity  for  each  kernel.  The 
derived  kinetic  energy  curves  show  a  ~ 30-fold  increase  in 


kinetic  energy  during  the  flare  with  a  decay  time  similar  to  the 
decay  rate  of  the  Ha  intensity.  The  measure  of  a  kernel’s  kinetic 
energy  is  imprecise  because  the  measured  motion  is  apparent 
motion  of  the  underlying  plasma.  In  the  model  of  two  ribbon 
solar  flares  put  forth  by  Priest  &  Forbes  (2002),  the  apparent 
velocity  is  a  better  indicator  of  rates  of  magnetic  reconnection 
rather  than  of  plasma  motion.  Despite  this  caveat,  the  derived 
kernel  kinetic  energy  is  a  useful  measure  because  it  combines  the 
size  of  the  kernels  with  their  apparent  velocities  to  characterize 
the  flare’s  evolution,  thus  representing  the  two  in  a  quantitative 
way. 

The  non-thermal  25-50  keV  emission  between  16:42  and 
16:43  UT  from  RHESSI  (Lin  et  al.  2002)  of  the  May  13  flare 
is  contoured  in  Figure  1 1  over  an  Ha  image  taken  at  16:42  UT. 
The  high-energy  emission  is  centered  over  the  flare  ribbons  but 
is  discontinuous  across  the  ribbon.  There  are  three  localized 
points  from  which  the  majority  of  the  X-ray  emission  comes. 
The  parts  of  the  flare  ribbons  exhibiting  the  most  displacement 
show  generally  lower  X-ray  intensity  than  the  more  stationary 
segments  of  the  flare  ribbons. 

4.  SCB  PROPERTIES 

The  measured  properties  of  SCBs  can  be  studied  using  two 
scenarios:  each  SCB  kernel  is  considered  as  independent  and 
isolated  compact  brightenings;  each  compact  brightening  is 
considered  to  be  fragments  of  a  dynamic  system  that  has  larger 
structure.  This  is  similar  to  the  approach  used  for  flares.  Again, 
considering  both  scenarios  has  its  benefits.  Identifying  kernels  as 
independent  elements  reveals  structural  differences  in  compact 
brightenings  and  is  discussed  in  Section  4.1.  Associating  kernels 
with  their  contextual  surroundings  provides  some  insight  into 
the  physical  structure  causing  SCBs.  An  aggregate  assessment  of 
SCBs  is  addressed  in  Section  4.2.  To  mitigate  the  effects  of  false 
or  marginal  detections,  only  SCBs  with  intensities  two  standard 
deviations  above  the  mean  background  intensity  are  considered 
for  characterizing  the  nature  of  the  ephemeral  brightening. 

4.1.  Qualities  of  Individual  SCBs 

SCBs,  although  related  to  erupting  flare  ribbons,  are  distinctly 
different  from  the  flare  kernels  discussed  in  Section  3.  Six  SCBs 
are  chosen  from  the  2005  May  13  event  as  an  example  of  these 
ephemeral  phenomena.  The  locations  of  these  six  events  are 
letters  A-F  in  Figure  6.  The  Ha  normalized  intensity  of  each 
of  these  kernels  is  shown  in  Figure  7.  The  SCB  intensity  curve 
is  significantly  different  from  the  flare  kernel  curves  shown  on 
the  left  side  of  the  figure.  SCB  curves  are  impulsive;  they  have 
a  sharp  peak  and  then  return  to  background  intensity  in  the 
span  of  about  12  minutes.  Nearly  all  of  the  SCBs  shown  here 
peak  in  intensity  before  the  peak  of  the  flare  intensity  curve, 
shown  as  a  vertical  dashed  line  in  Figure  7.  SCBs  B  and  E 
both  appear  to  have  more  internal  structure  than  the  other  SCBs 
and  last  noticeably  longer.  This  is  most  likely  caused  by  several 
unresolved  SCB  events  occurring  in  succession. 

Nearly  all  SCBs  have  Ho'  line  center  intensities  40%-60% 
above  the  mean  intensity  of  the  quiet  chromosphere  (Figure  12). 
From  this  increased  initial  intensity,  SCBs  brighten  to  intensities 
75%-130%  above  the  quiet  Sun. 

Examining  Doppler  velocity  measurements  from  the  SCB 
locations  reveals  three  distinct  types  of  SCBs  (Figure  12).  A 
type  la  SCB  has  an  impulsive  intensity  profile  and  an  impulsive 
negative  Doppler  profile  that  occurs  simultaneously  or  a  few 
minutes  after  the  peak  brightening.  In  this  study,  a  negative 
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Figure  10.  Integrated  unsigned  propagation  speed,  acceleration,  and  energy  of  flare  kernels  as  a  function  of  time  for  each  of  the  three  events  studied.  The  vertical 
dashed  line  indicates  the  maximum  flare  intensity  for  the  flaring  event.  The  first  column  shows  the  integrated  speed  curves.  The  second  column  shows  the  acceleration 
of  the  flare  ribbons.  The  third  column  shows  the  integrated  kinetic  energy.  Equation  (8),  as  a  function  of  time  for  each  of  the  three  events  studied.  The  2004  flare  had 
significant  noise  in  the  original  data  and  a  location  on  the  limb  contributing  to  the  inconsistent  measurements. 
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Figure  11.  Ho;  image  of  the  2005  May  13  flare  at  peak  intensity  is  shown  with 
the  reconstructed  RHESSI  25-50  keV  pass  band  integrated  between  16:42  and 
16:43  UT  contoured  on  top  at  10%,  30%,  50%,  70%,  and  90%  of  the  maximum 
intensity.  The  held  of  view  of  the  RHESSI  image  [—220:— 95,  180:305]  is 
confined  within  the  flaring  region. 


velocity  is  associated  with  motion  away  from  the  observer  and 
into  the  Sun.  A  type  Ib  SCB  has  a  similar  intensity  and  Doppler 
profile  as  a  type  la  but  the  timing  of  the  impulsive  negative 
velocity  occurs  several  minutes  before  the  peak  intensity  of 
the  SCB.  The  type  Ib  SCB  shown  in  Figure  12  has  a  negative 
velocity  that  peaks  10  minutes  before  the  peak  intensity  and 
returns  to  a  stationary  state  before  the  Ha  line  center  intensity 
has  decayed.  A  type  II  SCB  has  a  positive  Doppler  shift 
perturbation  that  often  lasts  longer  than  the  emission  in  the  Ha 
intensity  profile.  The  timing  of  both  are  nearly  coincidental.  A 
type  III  SCB  demonstrates  variable  dynamics.  It  has  a  broad  Ha 
intensity  line  center  with  significant  substructure.  A  type  III  SCB 
begins  with  a  negative  Doppler  profile  much  like  a  type  I.  Before 
the  negative  velocity  perturbation  can  decay  back  to  continuum 
levels,  there  is  a  dramatic  positive  velocity  shift  within  three 
minutes  with  an  associated  line  center  brightening.  In  all  types 
of  SCBs  the  typical  magnitude  of  Doppler  velocity  perturbation 
is  between  2  and  5  km  s“^  in  either  direction  perpendicular  to 
the  solar  surface. 

Of  all  off-flare  compact  brightenings  detected  using  the  au¬ 
tomated  techniques,  59%  of  SCBs  have  no  discernible  Doppler 
velocity  in  the  May  6  event  and  21%  in  the  May  13  event.  Out 
of  the  SCBs  that  do  have  an  associated  Doppler  velocity,  the 
May  6  event  has  35%  of  type  I,  54%  of  type  II,  and  11%  of 
type  III.  The  May  13  event  has  52%  of  type  I,  41%  of  type  II, 
and  17%  of  type  III.  When  totaling  all  of  the  Doppler  velocity 
associated  SCBs  between  both  events,  41%  are  of  type  I,  46% 
of  type  II,  and  13%  of  type  III.  Both  line  center  brightenings  and 


8 

Approved  for  public  release;  distribution  is  unlimited. 


The  Astrophysical  Journal,  750:145  (13pp),  2012  May  10 


Kirk  et  al. 


Figure  12.  Line  center  intensity  and  Doppler  velocity  measurements  for  three  different  types  of  SCBs  observed  in  both  2005  flares.  The  vertical  dashed  line  marks  the 
peak  Ho;  intensity  of  the  associated  SCB.  A  negative  Doppler  velocity  is  away  from  the  observer  and  into  the  Sun. 


SCB  Duration 


Figure  13.  Histogram  of  the  duration  of  all  SCB  detections  in  all  three  events. 
The  SCBs  have  a  measured  mean  duration  of  5.7  minutes  (dot-dashed  line)  and 
a  median  duration  of  3.1  minutes  (dashed  line). 


Doppler  velocities  are  filtered  such  that  positive  detections  have 
at  least  a  two  standard  deviation  peak  above  the  background 
noise  determined  in  the  candidate  detection  (Kirk  et  al.  201 1). 

4.2.  SCBs  in  Aggregate 

Examining  SCBs  as  a  total  population,  the  intensity  bright- 
enings  have  a  median  duration  of  3.1  minutes  and  a  mean  du¬ 
ration  is  5.7  minutes  (Figure  13).  The  duration  is  characterized 
by  the  FWHM  of  the  SCB  intensity  curve  (examples  of  these 
curves  are  shown  in  Figure  6).  A  histogram  of  the  distribution 
of  the  number  of  SCB  events  as  a  function  of  duration  shows  an 


exponential  decline  in  the  number  of  SCB  events  between  2  and 
30  minutes  (Figure  13).  The  duration  of  SCBs  is  uncorrelated 
with  both  distance  from  flare  center  and  the  peak  intensity  of 
the  SCB. 

Examining  individual  SCBs’  distance  from  the  flare  center 
as  a  function  of  time  provides  a  way  to  extract  the  propagation 
trends  of  SCBs  around  the  flare  (Figure  14).  The  bulk  SCBs 
appear  between  1.2  and  2.5  Mm  away  from  flare  center  in 
a  1  hr  time  window.  At  the  extremes,  SCBs  are  observed  at 
distances  up  to  ~ 5  Mm  and  several  hours  after  the  flare  intensity 
maximum.  On  the  top  two  panels  of  Figure  14,  the  shade  of  the 
mark  corresponds  to  the  intensity  of  the  SCB  measured  where 
the  lighter  the  mark,  the  brighter  the  SCB.  The  center  of  the  flare 
is  determined  by  first  retaining  pixels  in  the  normalized  data  set 
above  1.35.  Then  all  of  the  images  are  co-added  and  the  center 
of  mass  of  that  co-added  image  is  set  to  be  the  “flare  center.”  The 
SCBs  tend  to  clump  together  in  the  time-distance  plot  in  both 
the  May  6  and  May  13  events.  Generally,  the  brighter  SCBs  are 
physically  closer  to  the  flare  and  temporally  occur  closer  to  the 
flare  peak.  This  intensity  correlation  is  weak  and  qualitatively 
related  to  distance  rather  than  time  of  brightening.  Statistics 
from  the  2004  November  9  event  are  dominated  by  noise. 

Identifying  SCBs  with  negative  Doppler  shifts  within 
three  minutes  of  the  Ha  peak  intensity  limits  SCBs  to  type  I 
only.  With  this  method  of  filtering,  some  trends  become  more 
apparent  in  a  time-distance  plot  of  SCBs  (Figure  14,  bottom 
panels).  In  this  view  of  SCBs,  they  tend  to  cluster  together  in 
time  as  well  as  distance,  with  two  larger  groups  dominating  the 
plots. 

To  fit  a  slope  to  these  populations,  a  forward-fitting  tech¬ 
nique  is  employed  similar  to  a  linear  discriminate  analysis. 
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Figure  14.  Derived  statistics  for  all  SCBs  for  the  2005  May  6  and  2005  May  13  events.  All  panels  show  the  distance  at  which  the  SCB  occurs  vs.  time  from  the  Ho; 
flare  peak.  The  color  of  each  plotted  point  is  representative  of  its  relative  intensity;  the  lighter  points  are  higher  intensity  detections.  Bottom:  the  same  plots  as  the  top 
two  panels  except  with  a  Doppler  filter  applied — each  point  had  to  have  an  SCB  Type  I  Doppler  shift  associated  with  the  line  center  brightening.  The  dashed  lines  are 
fit  to  the  two  visually  distinct  groups  in  each  event. 


This  method  requires  the  user  to  identify  the  number  of 
groups  to  be  fitted  and  the  location  of  each  group.  The  fit¬ 
ting  routine  then  searches  all  linear  combinations  of  features 
for  the  next  “best  point”  to  minimize  the  to  d  regression 
fit  of  the  candidate  group.  Repeating  this  method  over  all 
the  points  in  the  set  produces  an  ordered  set  of  points  that 
when  fit,  have  an  increasing  y}  value.  A  threshold  is  taken 
where  the  derivative  of  the  y}  curve  increases  to  beyond  one 
standard  deviation,  which  has  the  effect  of  identifying  where 
the  x^  begins  to  increase  dramatically.  Running  this  routine 
several  times  minimizes  the  effect  of  the  user  and  provides 
an  estaminet  of  the  error  associated  with  the  fitted  line.  This 
method  has  two  caveats.  First,  this  fitting  method  relies  on 
the  detections  having  Poissonian  noise.  This  is  not  necessar¬ 
ily  correct  since  the  detection  process  of  compact  brightenings 
introduces  a  selection  bias.  Second,  the  fitting  method  makes 
the  assumption  that  no  acceleration  occurs  in  the  propaga¬ 
tion  of  SCBs.  This  is  a  reasonable  approximation  but  from 
studies  of  Moreton  waves  in  the  chromosphere  (e.g.,  Bala- 
subramaniam  et  al.  2010),  a  constant  velocity  propagation  is 
unlikely. 

Applying  the  forward-fitting  technique  these  two  data  sets 
yields  two  propagation  speeds:  a  fast  and  a  slow  group.  The 
May  6  event  has  propagation  speeds  of  71  zb  2  km  s“^  and 
49  zb  2  km  s“^  The  May  13  event  has  propagation  speeds  of 
106  zb  7  km  s“^  and  33  zb  4  km  s“^ .  Beyond  these  detections,  the 
May  1 3  event  has  a  third  ambiguous  detection  with  a  propagation 


speed  of  —7  zb  5  km  s“^ .  The  implications  of  these  statistics  are 
discussed  in  Section  6. 

The  X-ray  intensity  measured  with  RHESSI  temporally  covers 
a  part  of  the  flaring  event  on  2005  May  1 3 .  RHESSI  data  coverage 
begins  at  16:37  UT  during  the  impulsive  phase  of  the  flare. 
Comparing  the  integrated  X-ray  intensity  curve  between  25 
and  50  keV  to  the  aggregate  of  SCB  intensities  integrated  over 
each  minute  yields  some  similarities  (Figure  15).  The  X-ray 
intensity  peaks  about  a  minute  after  the  integrated  SCBs  reach 
their  maximum  intensity.  The  decay  of  both  the  X-ray  intensity 
and  SCB  integrated  curve  occurs  approximately  on  the  same 
timescale  of  ~50  minutes. 

5.  PHYSICAL  INTERPRETATION  OF  RESULTS 

From  the  results  of  the  tracking  method,  we  can  gain  insight 
into  the  physical  processes  associated  with  SCBs  and  flares. 
First,  the  flare  ribbons  studied  are  both  spatially  and  temporally 
clumped  into  discrete  kernels.  The  relatively  smooth  ribbon 
motions  observed  in  Ho'  and  the  nearly  ideal  exponential  decay 
characterized  in  X-ray  flux  measurements  belie  the  substructure 
that  makes  up  a  two-ribbon  flare.  This  is  discussed  in  Section  5.1. 
Second,  the  brightening  in  the  flare  ribbons  is  caused  by 
a  distinctly  different  physical  process  than  the  one  causing 
SCBs.  In  this  sense,  SCBs  are  different  than  micro-flares.  Our 
interpretation  of  the  physical  origin  of  SCBs  is  addressed  in 
Section  5.2. 
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Integrated  SCB  Kernel  Intesity 


Figure  15.  Plot  showing  the  RHESSI  non-thermal  25-50  keV  intensity  in  black  and  the  Ho;  SCB  integrated  intensity  curve  in  gray  for  the  2005  May  13  flare.  The 
vertical  dashed  line  indicates  the  peak  of  the  integrated  SCB  curve. 


5.1.  Two  Ribbon  Chromospheric  Flares 

Flare  kernels  are  observed  to  appear  and  disappear  as  the 
underlying  flare  ribbons  evolve,  as  discussed  in  Section  3.1. 
Examining  individual  kernel  structure  suggests  that  there  is 
substructure  within  a  flare  ribbon  whose  elements  impulsively 
brighten  and  dim  within  the  brightness  that  encompasses  the 
intensity  curve.  These  results  support  the  premise  that  flares 
are  made  up  of  several  magnetic  fleld  lines  reconnecting  (e.g., 
Priest  &  Forbes  2002).  There  is  no  evidence  to  claim  that  an 
individual  flare  kernel  is  directly  tracking  a  loop  footpoint. 
Within  a  tracked  flare  kernel,  multiple  coronal  reconnection 
events  are  superimposed  to  produce  the  observed  asymmetries 
a  single  flare  kernel’s  light  curve. 

Maury  a  &  Ambastha  (2010)  tracked  subsections  of  the  2003 
October  28  X17  flare  as  it  evolved.  They  reported  peak  speeds 
ranging  from  ~10to60  km  s“^  (depending  on  the  spatially 
tracked  part  of  the  flare)  over  the  observed  span  of  13  minutes. 
Maurya  &  Ambastha  (2010)  also  reported  that  the  total  apparent 
distances  the  ribbons  traveled  were  "^lO'^  km.  The  October 
28  flare  is  several  orders  of  magnitude  greater  in  GOES 
X-ray  intensity  than  the  flares  considered  for  this  study.  In 
the  present  study,  despite  this  difference,  peak  speeds  of  flare 
kernels  observed  in  this  work  are  measured  at  ~2.3  km  s“^ 
and  the  mean  velocity  of  all  flare  kernels  is  ~0.2  km  s“^.  The 
maximum  distance  flare  kernels  traversed  was  ~2  x  10^  km  and 
the  average  distance  traveled  was  x  10^  km;  both  of  these 
values  are  similar  to  the  2003  October  28  event.  Although  the 
velocities  of  the  flare  ribbons  studied  here  are  at  least  an  order  of 
magnitude  less,  the  two  to  three  orders-of-magnitude  difference 
in  the  GOES  peak  intensity  between  the  October  28  flare  and 
the  flares  studied  here  implies  that  the  apparent  velocities  and 
distance  traveled  by  flare  ribbons  do  not  scale  linearly  with  the 
strength  of  the  eruption. 

The  bulk  of  the  apparent  motion  and  acceleration  in  the 
flare  ribbons  are  observed  after  the  peak  of  the  flare  intensity. 
Integrating  this  into  the  dynamical  model  of  a  two-ribbon  flare 
(Demoulin  &  Priest  1988)  implies  that  the  peak  intensity  of  the 
flare  occurs  in  the  low-lying  arcade  and  loses  intensity  as  the 
x-point  reconnection  progresses  vertically  to  higher  levels 
(Baker  et  al.  2009).  The  measured  divergence  of  the  ribbons 
dominates  the  motion  but  there  is  signiflcant  apparent  flow 
tangential  to  the  flare  ribbons.  This  predicates  that  there  is  lateral 
propagation  to  the  x-point  as  well  as  vertical  propagation. 


5.2.  A  Physical  Model  of  SCB s 

Figure  16  is  a  proposed  phenomenological  model  of  the 
overlying  physical  topology.  SCBs  are  hypothesized  to  be 
caused  by  electron  beam  heating  conflned  by  magnetic  loop 
lines  overarching  flare  ribbons.  These  overarching  loops  are 
analogous  to  the  higher  lying,  unsheared  tethers  in  the  breakout 
model  of  CMEs  (Antiochos  et  al.  1999).  As  the  flare  erupts, 
magnetic  reconnection  begins  in  a  coronal  x-point.  The  CME 
escapes  into  interplanetary  space,  the  remaining  loop  arcade 
produces  a  two-ribbon  flare,  and  the  tethers  reconnect  to  a 
new  equilibrium  position.  The  tether  reconnection  accelerates 
trapped  plasma  which  impacts  the  denser  chromosphere  causing 
observed  brightening.  This  description  implies  that  the  driver  of 
SCBs  is  the  eruption  of  an  CME.  Balasubramaniam  et  al.  (2006) 
found  a  strong  correlation  between  CMEs  and  the  presence  of 
SCBs. 

A  simple  loop  conflguration  without  localized  diffusion  or 
anomalous  resistance  implies  that  the  length  of  the  loop  directly 
proportional  to  the  travel  time  of  the  electron  beam.  This  means 
that  the  different  propagation  groups  observed  in  Figure  14 
result  from  different  physical  orientations  of  overarching  loops. 
The  time  it  takes  to  observe  a  chromospheric  brightening  can 
therefore  be  described  as 

L 

"^^scB  =  K  — ,  (9) 

Ve 

where  tscb  is  the  time  it  takes  for  accelerated  plasma  to  travel 
from  the  coronal  x-point,  along  the  magnetic  loop  lines,  impact 
the  chromosphere,  and  produce  a  brightening,  L  is  the  length 
of  the  tether  before  reconnection,  Ve  is  the  electron  velocity 
along  the  loop  line,  and  y  is  a  function  of  changing  diffusion 
and  resistance  between  chromosphere  and  corona.  In  physical 
scenarios,  Ve  is  likely  a  combination  of  the  Alfven  speed  and 
the  electron  thermal  velocity. 

The  slopes  plotted  in  Figure  14  show  two  different  propa¬ 
gation  speeds  for  SCBs:  33-49  km  s“^  and  71-106  km  s“^. 
The  sound  speed  in  the  chromosphere  is  approximately 
Q  ~  10  km  s“^  (Nagashima  et  al.  2009),  while  the  Alfven 
speed  in  the  upper  chromosphere  is  approximated  to  be  be¬ 
tween  ua  ~  10  and  100  km  s~^  (Aschwanden  2005).  Both  of 
the  SCB  propagation  speeds  are  signiflcantly  above  the  approx¬ 
imated  sound  speed;  however,  they  fall  reasonably  well  into  the 
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Figure  16.  Schematic  model  of  the  physical  cause  of  SCBs.  A  dashed  line  marks  the  neutral  line  in  this  model.  The  central  loops  represent  the  emerging  two-ribbon 
flare  arcade.  The  arches  represent  held  lines  connecting  outside  of  the  flare  ribbons.  The  arrows  on  top  of  both  the  loops  and  arches  show  the  direction  trapped  plasma 
flows  as  the  flare  begins  to  erupt  and  the  stars  suggest  the  location  where  the  plasma  is  disrupted.  Once  perturbed,  the  trapped  plasma  streams  down  the  loop  lines 
and  impacts  the  chromosphere  causing  the  Ho;  intensity  brightening.  The  different  orientations  of  off-flare  loops  account  for  the  different  propagation  speeds  of  SCBs 
observed. 
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Figure  17.  Diagram  of  the  physical  dynamics  occurring  in  a  single  SCB.  Electrons  and  protons  accelerated  by  magnetic  reconnection  further  up  in  the  corona  come 
streaming  down  the  flare  loop.  Since  the  mean  free  path  of  electrons  is  signiflcantly  less,  the  electrons  deposit  energy  into  the  upper  chromosphere.  Unable  to  radiate 
the  energy  as  heat,  the  chromosphere  explosively  responds  by  sending  material  back  up  the  flux  tube  (chromospheric  evaporation)  and  a  reaction  shock  propagates 
toward  the  photosphere.  This  reaction  shock  explains  the  velocities  observed  in  Type  I  SCBs. 


range  of  the  Alfven  speed.  Assuming  the  propagation  speed 
of  non-thermal  plasma  along  the  loop  line  is  consistent  be¬ 
tween  flare  loops,  the  different  propagation  speeds  would  then 
imply  different  populations  of  loops  being  heated  as  the  flare 
erupts. 

Examining  single  SCBs,  the  coincident  Doppler  recoil  with 
the  Hckf  intensity  presents  a  contradiction.  If  an  SCB  is  an 
example  of  compact  chromospheric  evaporation  (e.g.,  Dennis 
&  Schwartz  1989),  then  the  only  Doppler  motion  should  be 
outward,  the  opposite  of  what  is  observed.  Figure  17  presents 
a  possible  solution  to  this.  Since  the  scattering  length  of 
electrons  is  significantly  smaller  than  that  of  protons,  the 


electron  beam  impacts  the  mid-chromosphere  and  deposits 
energy  into  the  surrounding  plasma  while  the  protons  penetrate 
deeper.  This  deposited  heat  cannot  dissipate  effectively  through 
conduction  or  radiation  and  thus  expands  upward  into  the  flux 
tube.  To  achieve  an  expansion  with  rates  of  a  few  km  s“^ 
as  observed,  the  chromospheric  heating  rate  must  be  below 
Eh  ^  10^^  erg  cm  ^  s  ^  (Fisher  et  al.  1984).  As  a  reaction 
to  this  expansion,  a  reaction  wave  propagates  toward  the  solar 
surface  in  the  opposite  direction  of  the  ejected  plasma.  Since  the 
Doppler  measurements  are  made  in  the  wings  of  the  Ha  line, 
the  location  of  the  Doppler  measurement  is  physically  closer  to 
the  photosphere  than  the  Ha  line  center.  Thus,  the  observer  sees 
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the  heating  of  the  Ha  line  center  and  coincidentally  observes 
the  recoil  in  the  lower  chromosphere. 

In  type  II  SCBs  an  upflow  is  observed.  This  is  an  example  of 
a  classic  chromospheric  evaporation  where  the  heated  plasma 
in  the  bulk  of  the  chromosphere  is  heated  and  ablated  back 
up  the  flux  tube.  In  contrast,  the  type  III  class  of  SCB  shown 
in  Figure  12  presents  an  interesting  anomaly  to  both  the  other 
observed  SCBs  and  the  model  proposed  in  Figure  17.  Since 
there  is  an  initial  down-flow,  the  beginning  state  of  type  III  SCBs 
are  similar  to  type  Is.  As  the  recoil  is  propagated,  a  continual 
bombardment  of  excited  plasma  (both  protons  and  electrons) 
impacts  the  lower  chromosphere,  causing  ablation  and  changing 
the  direction  of  flow. 

6.  SUMMARY  AND  CONCLUSIONS 

Two  conclusions  can  be  made  about  the  flares  studied  in  this 
work:  first,  the  asymmetrical  motions  of  the  flare  ribbons  imply 
that  the  peak  flare  energy  occurs  in  the  low-lying  arcade,  and 
second,  flare-related  SCBs  appear  at  distances  on  the  order  of 
10^  km  and  have  properties  of  sites  of  compact  chromospheric 
evaporation.  Maury  a  &  Ambastha  (2010)  estimated  reconnec¬ 
tion  rates  from  a  measurement  of  the  photospheric  magnetic 
field  and  the  apparent  velocity  of  chromospheric  flare  ribbons. 
Therefore,  it  would  be  possible  to  estimate  the  reconnection  rates 
using  this  technique  with  the  addition  of  photospheric  magne¬ 
tograms.  Associating  vector  magnetograms  with  this  technique 
and  a  careful  consideration  of  the  observed  Doppler  motions 
underneath  the  ribbons  would  provide  a  full  three-dimensional 
method  for  estimating  the  Lorentz  force  for  subsections  of  a 
flare  ribbon. 

SCBs  originate  during  the  impulsive  rise  phase  of  the  flare  and 
often  precede  the  Ha  flare  peak.  SCBs  are  found  to  appear  with 
similar  qualities  as  compact  chromospheric  ablation  confirming 
the  results  of  Pevtsov  et  al.  (2007).  The  heuristic  model  pro¬ 
posed  in  Figure  17  requires  unipolar  magnetic  field  underneath 
SCBs.  The  integration  of  high-resolution  magnetograms  and  a 
field  extrapolation  into  the  chromosphere  would  confirm  this 
chromospheric  evaporation  model.  One  consequence  of  chro¬ 
mospheric  evaporation  is  that  the  brightening  should  also  be 
visible  in  EUV  and  X-ray  observations  due  to  non-thermal  pro¬ 
tons  interacting  with  the  lower  chromosphere  or  photosphere. 

SCBs  are  a  special  case  of  chromospheric  compact  brighten¬ 
ing  that  occur  in  conjunction  with  flares.  The  distinct  nature  of 
SCBs  arises  from  their  impulsive  brightenings,  unique  Doppler 
velocity  profiles,  and  origin  in  the  impulsive  phase  of  flare  erup¬ 
tion.  These  facts  combined  demonstrate  that  SCBs  have  a  non- 
localized  area  of  influence  and  are  indicative  of  the  conditions 


of  the  entire  flaring  region.  They  can  possibly  be  understood 
by  a  mechanism  in  which  a  destabilized  overlying  magnetic 
arcade  accelerates  electrons  along  magnetic  tubes  that  impact 
a  denser  chromosphere  to  result  in  an  SCB.  This  distinguishes 
SCBs  from  the  flare  with  which  they  are  associated.  In  a  future 
work,  we  plan  to  address  physical  mechanisms  to  explain  the 
energetic  differences  between  SCBs  and  flares  to  explain  the 
coupled  phenomena. 

The  authors  thank  USAF/AFRL  grant  FA9453- 11-1-0259; 
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